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1. INTRODUCTION 

In recent years nucleon properties have been in the focus of manifestly Lorentz covariant 
Chiral Perturbation Theory (ChPT), improved lattice QCD computations and chiral 
extrapolations (see e.g. Refs. p]-[23]). The lattice formulation of QCD is well established 
and is a powerful tool for studying the structure of nucleons. The computation of 
nucleon properties in lattice QCD is progressing with steadily increasing accuracy [I]- 
[H]. Accurate computations of the nucleon mass with dynamical fermions and two 
active flavors are now possible [9j [10] in lattice QCD. In practice, these computations 
are so far limited to relatively large quark masses. Direct simulations of QCD for 
light current quark masses, near the chiral limit, remain computationally intensive. 
To extract predictions for observables, lattice data generated at high current quark 
masses have to be extrapolated to the point of physical quark or pion mass. Therefore, 
one of the current aims in lattice QCD is to establish the quark mass dependence of 
quantities of physical interest, such as the nucleon mass, magnetic moments and form 
factors. The major tool in establishing the current quark mass dependence of lattice 
QCD results are methods based on chiral effective field theory. Recent extrapolation 
studies of lattice results concern the nucleon mass [U El EED, (El [Hj, its axial vector 
coupling constant and magnetic moments [H]-[T3], the pion-nucleon sigma term, charge 
radii [15], form factors [16]- [H], and moments of structure functions [19] . The chiral 
expansion in chiral effective field theory (xEFT) has been used to study the quark 
mass (pion mass) dependence of the magnetic moments, magnetic form factors and the 
axial-vector coupling constant [201 2 1 J of the nucleon for extrapolations of lattice QCD 
results, so far determined at relatively large quark masses corresponding to pion masses 
of m n > 0.6 GeV, down to physical values of m n . In the chiral limit, with m, — > 0, QCD 
at low energies is realized in the form of an effective field theory with spontaneously 
broken chiral symmetry, with massless pions as the primary active degrees of freedom. 
The coupling of the chiral Goldstone bosons to these spin- 1/2 matter field produces the 
so-called "pion-cloud" of the nucleon, an important component of nucleon structure at 
low energy and low momentum scales. 

In the present paper we investigate the dependence of nucleon properties (mass, 
magnetic moments and electromagnetic form factors) on pseudoscalar meson masses 
applying the perturbative chiral quark model (PCQM) [21]- [26]. In the PCQM baryons 
are described by three relativistic valence quarks confined in a static potential, which 
are supplemented by a cloud of pseudoscalar Goldstone bosons, as required by chiral 
symmetry. This simple phenomenological model has already been successfully applied 
to the charge and magnetic form factors of baryons, sigma terms, ground state masses of 
baryons, the electromagnetic N —>■ A transition, and other baryon properties |24|-[26|. 
Note that in Refs. [271 EE] we extend this approach by constructing a framework which 
is manifestly Lorentz covariant and aims for consistency with ChPT. 

In this work, our strategy is as follows. First, we discuss the nucleon properties 
(mass, magnetic moments and electromagnetic form factors) in dependence on the 
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pion mass in the two-flavor sector. Second, we extend our formalism to the three- 
flavor sector including kaon and ^-meson degrees of freedom with fixed masses. All 
calculations are performed at one loop. The chiral limit, where current quark masses 
approach zero with rh, m s — > 0, is well defined. We compare the obtained quark mass 
dependence of the nucleon observables to the results of other approaches (lattice QCD, 
chiral extrapolations). 

The paper is organized as follows. In Sec. II we give a short overview of our 
approach. In Sec. Ill we the discuss dependence of nucleon properties on the variation 
of the pion mass in the two- and three-flavor picture in the context of the PCQM and 
compare them to other theoretical approaches. In Sec. IV we give our conclusions. 



2. THE PERTURB ATIVE CHIRAL QUARK MODEL 



The perturbative chiral quark model [23]- [26] is based on an effective chiral Lagrangian 
describing baryons by a core of three valence quarks, moving in a central Dirac field with 
V e f[(r) = S(r) + 7°V A (r), where r —\ x |. In order to respect chiral symmetry, a cloud 
of Goldstone bosons (n, K and n) is included, which are treated as small fluctuations 
around the three-quark core. The model Lagrangian is 

C(x) = ${x) [i $ - 7°V(r) - M] ip(x) 
F 2 r 

+ — Tr^ l U(x)d'"U\x) + 2MB{U{x) + U ] {x)) 

- i){x) S(r) I 2 — + 1 2 —\ip(x), (1) 

where ip = (u,d,s) is the triplet of quark fields, U = exp[z$/F] is the chiral field 
in the exponentional parametrization, F = 88 MeV is the pion decay constant in the 
chiral limit [29], A4 = diag{m u , m<2, m s } is the mass matrix of current quarks and 
B = — (0\uu\0)/F 2 = — (0\dd\0)/F 2 is the quark condensate constant. In the numerical 
calculations we restrict to the isospin symmetry limit m u = m<i = rh. We rely on the 
standard picture of chiral symmetry breaking [30] and for the masses of pseudoscalar 
mesons we use the leading term in their chiral expansion (i.e. linear in the current 
quark mass). By construction, our effective chiral Lagrangian is consistent with the 
known low-energy theorems (Gell-Mann-Okubo and Gell-Mann-Oakes-Renner relations, 
partial conservation of axial current (PCAC), Feynman-Hellmann relation between pion- 
nucleon a-term and the derivative of the nucleon mass, etc.). The electromagnetic field 
is included into the effective Lagrangian (0Q) using the standard procedure, i.e. the 
interaction of quarks and charged mesons with photons is introduced using minimal 
substitution. 

To derive the properties of baryons, which are modeled as bound states of valence 
quarks surrounded by a meson cloud, we formulate perturbation theory and restrict the 
quark states to the ground-state contribution with i[)(x) = boUo(x) exp(— i£ot), where 6 
is the corresponding single-quark annihilation operator. The quark wave function Uq(x) 
belongs to the basis of potential eigenstates used for expanding the quark field operator 
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ip{x). In our calculation of matrix elements, we project quark diagrams on the respective 
baryon states. The baryon states are conventionally set up by the product of the SU(6) 
spin-flavor and SU(3) C color wave functions, where the nonrelativistic single quark spin 
wave function is simply replaced by the relativistic solution Uq(x) of the Dirac equation 

-z 7 °7 ■ V + 7 °S(r) + V(r) - S ] u (x) = 0, (2) 

where So is the single-quark ground-state energy. 

For the description of baryon properties, we use the effective potential V e s{r) with 
a quadratic radial dependence [Ml EH] : 

S{r) =M 1 + Cl r 2 , V(r) = M 2 + c 2 r 2 (3) 

with the particular choice 

1 - 3p 2 1 + 3p 2 p 

Here, R and p are parameters related to the ground-state quark wave function u : 

1 



Uq(x; i) = N exp 



x 2 



2R 2 



Xs(i)Xf(i)Xc(i) , (5) 



ipa(i) ■ x/R 

where JV = [vr 3 / 2 /^! + 3p 2 /2)]- 1 / 2 is a normalization constant; Xs, Xfi Xc are the 
spin, flavor and color quark wave functions, respectively. The index "z" stands for the 
i-th quark. The constant part of the scalar potential Mi can be interpreted as the 
constituent mass of the quark, which is simply the displacement of the current quark 
mass due to the potential S(r). The parameter p is related to the axial charge qa of the 
nucleon calculated in zeroth-order (or 3q-core) approximation: 



3 V 1 + § P 2 

Therefore, p can be replaced by qa using the matching condition (6). The parameter R 
is related to the charge radius of the proton in the zeroth-order approximation as 

(rl)lo = j d?xul{x)x 2 uo{x) = 3J f 1 j^ 2 - (7) 

In our calculations we use the value ^=1.25. Therefore, we have only one free parameter 
in our model, that is R or (rj^^Q. In previous publications R was varied in the region 
from 0.55 fm to 0.65 fm, which corresponds to a change of (r|;}£ from 0.5 to 0.7 fm 2 . 
Note that for the given form of the effective potential ([3]) the Dirac equation ([2]) can be 
solved analytically [for the ground state see Eq.fjSJ), for excited states see Ref. [26]]. 
The expectation value of an operator A is then set up as: 

(A)= B (<Po\f2^ [ d 4 Xl ... f d 4 x n T[£ I {x 1 )...£ I {x n )A)\cf> )*, (8) 

n=l n ' J J 

where the state vector \<f) > corresponds to the unperturbed three-quark state (3g- 
core). Superscript "B" in the equation indicates that the matrix elements have to be 
projected onto the respective baryon states, whereas subscript "c" refers to contributions 
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from connected graphs only. Here Cj{x) is the appropriate interaction Lagrangian. 
For the purpose of the present paper, we include in Ci(x) the linearized coupling of 
pseudoscalar fields with quarks and the corresponding coupling of quarks and mesons 
to the electromagnetic field (see details in Ref. [21]): 

d(x) = - ^{x)ij°^§^S{r)tfj{x) - eA^{x)tfj{x)YQ^{x) 
b 

8 , 

- eA,(x) J2 [f*S + ;§] *i(x)P* s (x) + ■■■, (9) 

i,j=l V 

where fijk are the SU(3) antisymmetric structure constants. 

For the evaluation of Eq. ([8]) we apply Wick's theorem with the appropriate 
propagators for the quarks and pions. For the quark propagator we use the vacuum 
Feynman propagator for a fermion in a binding potential restricted to the ground-state 
quark wave function with 

iG (x, y) = u (x) u (y) e ^o-y ) g ^ Q _ yo) (10) 
For the meson field we use the free Feynman propagator for a boson field with 

iA^x -y) = (0\T{^(x)^(y)}\0) = [ e~ ik ^ A*(fc) , (11) 



(27T) 

where A*(fc) = [M| — k 2 — i0 + ] _1 is the meson propagator in momentum space and M$ 
is the meson mass. 

The physical nucleon mass at one loop is given by 

m N = m c ° re + Am N (12) 

where 

m c ° re = 3{S + 7 m} = 3^ + ^M 2 } (13) 

is the contribution of the three-quark core (the second term in the r.h.s. of Eq. f|T3|) is 
the contribution of the current quark mass) and 

Am N = IT MC + IT ME (14) 

is the nucleon mass shift due to the meson cloud contribution. The diagrams that 
contribute to the nucleon mass shift Am^ at one loop are shown in Fig.l (see details 
in Refs. [24~| 125]). Fig.la corresponds to the so-called meson-cloud (MC) contribution 
and Fig. lb is the meson-exchange (ME) contribution. The operators n MC and n ME are 
functions of the meson masses and are expressed in terms of the universal self-energy 
operator 

n(M|) = -J| 4 , ^ = 7T,K,rj. (15) 

Here we introduce a notation for the structure integral in terms of which all further 
formulas can be expressed: 

oo 

ir = (^) 2 / dP ^ NN S\ M,N = 0,1,2,... (16) 
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jQN — T — 



dp P N F 2 NN ( P 2 



x r 



9a 

7lF 

N+l\/l N + l 



9 r> N+l 



(3 + 



(N+l)(N + 3) 



17) 



.2 4' 32 
where (3 = 2p 2 / (2 — p 2 ) . The function F n NN(p 2 ) is the nNN form factor normalized to 
unity at zero recoil (p 2 = 0): 

Q 2 R 2 \ r g 2 i? 2 



irN N 



cxp 



4 



1 



The meson cloud contributions to the mass shift are then given 
1) in SU(2) as 

\l 



■qMC . 
■qME 

2) in SU(3) as 

j-jMC . 
■qME . 



400 
90 

400 
81 



U(M- 



7r) ' 



U(M- 



7r) > 



(19) 



54 



n(M 2 ) + Ti(Ml) 

400 V *' 400 V K; 



9 
4~00 



n ( M ')> 



90 



6 



— n(M 2 ) n(M 2 ) 

400 V 400 V ^ 



(20) 



Exact expressions for the nucleon electromagnetic form factors can be found in Ref. |24J. 
Here we just present the typical results for the magnetic moments in SU(2) and SU(3). 
The two-flavor result is obtained from the three-flavor one when neglecting kaon and 
rj- meson contributions. The magnetic moments of the nucleons, p p and p n , are given by 
the expressions 



s 



--u LO 



1 + 5 



400 



21/f + 21 J* + 4/f 



50 



K f ) 



where 



LO _ 2m N pR 



1 + fp 2 



is the leading-order contribution to the proton magnetic moment. The factor 

n MC i + |p 2 \ 2 



m + 



1-|P 2 



Ip 2 



i + Ip 2 



(21) 
(22) 

(23) 



defines the NLO correction to the nucleon magnetic moments due to the modification 
of the quark wave function |24j). 
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3. NUMERICAL RESULTS 

In this section we discuss the numerical results for the dependence of nucleon properties 
on a variation of the pseudoscalar meson masses. 

In Table 1 we present our results for the nucleon mass in dependence on the 
pion mass. Both the SU(2) version, considering only the pion cloud contribution, and 
the SU(3) variant, including in addition kaon and 77-meson cloud contributions, are 
indicated. The total result is normalized to the physical value (coinciding with the 
proton mass treated as the reference point) = m p = 938.27 MeV by fixing the 
ground-state quark energy to So ~ 397 MeV [in case of SU(2)] and So ~ 411 MeV [in 
case of SU(3)]. We also indicate the separate contributions of the 3q-core and the meson 
cloud, and, in addition, the value obtained in the chiral limit, consistent with the values 
of Refs. [271 I3T1 [14"] . For the dependence on the pion mass we choose mass values in the 
range of M 2 ~ 0.15 — 1.2 GeV 2 and the resulting nucleon mass is directly compared to 
either chiral extrapolations or lattice data [3, E]- In both comparisons our results are 
consistent with the corresponding values of either the extrapolations or the lattice data. 
In Figs. 3 and 4 we indicate the full functional dependence of the nucleon mass on M 2 , 
M\ and M 2 , respectively, and compare them (in case of the SU(2) M^-dependence) 
to results of lattice QCD at order p 3 and p A of Ref. [I]. In Fig. 5 we compare the 
results for the nucleon mass both in SU(2) and SU(3) to lattice QCD data from various 
collaborations [H [8] as functions of M 2 . 

In Table 2 we present our results for the nucleon magnetic moments in the SU(2) 
version for different values of the model scale parameter R at M v = and at the physical 
pion mass M^ hys . In Table 3 we give the analogous results for SU(3). In Figs. 6 and 7 
we draw the curves for the nucleon magnetic moments as functions of M 2 and compare 
them to results of lattice QCD [131 EH]- Note, that the nucleon magnetic moments are 
not sensitive to a variation of the strange current quark mass m s . In Figs. 8 and 9 we 
demonstrate the sensitivity of the nucleon magnetic moments as functions of M 2 on the 
variation of the scale parameter R. 

Finally, in Figs. 10-12 we present results for the nucleon form factors G P E (Q 2 ), 
G P M {Q 2 ) I Hp and G^iQ 2 ) / \i n at different values of M 2 . A larger pion mass leads to 
an increase of the normalized (at Q 2 = 0) form factors. 

4. CONCLUSIONS 

In this work we apply the perturbative chiral quark model at one loop to describe 
the dependence of nucleon properties on the meson masses. It has been previously 
verified that the model is successful in the explanation of many aspects of nucleon 
properties [23-25] , such as magnetic moments, the axial vector form factor, the N — ► A 
transition amplitude, the meson-nucleon sigma-term and ttN nucleon scattering. Here 
we demonstrate that the meson mass dependence of nucleon properties such as the 
mass, magnetic moments, electromagnetic form factors, both for the two and three flavor 
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variants, is reasonably described in comparison to present lattice data and extrapolations 
of these results. Given also the simplicity of this model approach, the evaluation at one 
loop seems sufficient to correctly describe the pion mass dependence of the discussed 
observables. 
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Table 1. Nucleon mass (in GeV) in the SU(2) and SU(3) versions and at different 
values for M% 
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Table 2. Nucleon magnetic moments in SU(2) 
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Table 3. Nucleon magnetic moments in SU(3) for fixed masses M\ and M% 
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Figure 1. Diagrams contributing to the nucleon mass shift. 
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Figure 2. Diagrams contributing to the nucleon electromagnetic form factors. 
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Figure 3. Dependence of nucleon mass wlm on meson mass M\: mjy(M^) in SU(2) 
(dotted line) and m N (M%), m N (M] ( ) m N (M?) in SU(3) (the other lines). 
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Figure 4. Dependence of nucleon mass tojv on pion mass M%: mjv(-M^) in SU(2) 
(dotted line), SU(3) (solid line) and from lattice QCD [4] (the others) at order p 3 
and v A - 



Dependence of nucleon properties on pseudoscalar meson masses 



14 




Figure 5. Nucleon dependencies mjy(M%) in SU(2) (dotted line) and in SU(3) (solid 
line) are compared to data from various collaborations as a function of M%. 
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Figure 6. M^-dependence of the proton magnetic moment /j, p (M%) to one-loop from 
sum rules [12] (dotted curve), lattice QCD [TH] (dashed curve) and our results in SU(3) 
(solid curve). 
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Figure 7. Af^-dependence of the neutron magnetic moment /x ra (Mj) to one-loop from 
sum rules [H] (dotted curve), lattice QCD [18] (dashed curve) and our results in SU(3) 
(solid curve). 
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(a) (b) 
Figure 8. fiN at M% = 0—1 GeV 2 and R = 0.5 - 0.7 fm in SU(2). 
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Figure 9. jj, N at M 2 = - 1 GeV 2 and R = 0.5 - 0.7 fm in SU(3). 
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Figure 10. Proton charge form factor G P E (Q 2 ) as function of M% in SU(2). 
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Figure 11. Proton magnetic form factor G P M (Q 2 ) as function of M% in SU(2). 
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Figure 12. Neutron magnetic form factor G^iQ 2 ) as function of M% in SU(2). 



